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Abstract. This paper describes a method for deriving approximate equations for irrota- 
tional water waves. The method is based on a 'relaxed' variational principle, i.e., on a 
Lagrangian involving as many variables as possible. This formulation is particularly suitable 
for the construction of approximate water wave models, since it allows more freedom while 
preserving a variational structure. The advantages of this relaxed formulation are illustrated 
with various examples in shallow and deep waters, as well as arbitrary depths. Using subordi- 
nate constraints (e.g., irrotationality or free surface impermeability) in various combinations, 
several model equations are derived, some being well-known, other being new. The models 
obtained are studied analytically and exact traveling wave solutions are constructed when 
possible. 



1. Introduction 

The water wave problem in fluid mechanics has been known since more than two hundreds 
years [13]. The classical mathematical formulation of surface gravity waves involves five 
equations: the irrotationality of the fluid flow, the fluid incompressibility, the bottom and 
the surface impermeabilities, and the surface isobarity [40] . This system of equations cannot 
generally be solved exactly and, historically, the water wave theory has been developed by 
constructing various approximations. In shallow water, we have the equations of Korteweg 
and de Vries [30], Boussinesq [7], Benjamin et al. [4], Serre [49], Green and Naghdi [23], 
Camassa and Holm [9], Degasperis-Procesi [14], and many other model equations. On finite 
depth and deep water, there is the celebrated nonlinear Schrodinger equation [40] and the 
equations of Dysthe [17], Trulsen et al. [55], Kraenkel et al. [31], among others. These 
equations are most often derived via some perturbation techniques and are thus valid for 
waves of small amplitude. Moreover, these equations are generally valid for a very limited 
range of the ratio wavelength / water depth and for narrow-banded spectra. However, for many 
applications it is necessary to use models uniformly valid for all depths and which are accurate 
for large amplitudes. It is well-known in theoretical physics that variational formulations are 
tools of choice to derive such approximations when small parameter expansions are inefficient. 

There are mainly two variational formulations for irrotational surface waves that are com- 
monly used, namely the Lagrangian of Luke [36] and the Hamiltonian of Petrov-Zakharov 
[46, 60]. Details on the variational formulations for surface waves can be found in review 
papers, e.g., [47, 48, 61]. The water wave problem is also known to have the multi-symplectic 
structure [8]. These variational principles have been exploited, in different variants, to build 
analytical and numerical approximations, e.g., [3, 28] just to mention a few references. 
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Luke's Lagrangian assumes that the flow is exactly irrotational, i.e., the Lagrangian in- 
volves a velocity potential but not explicitly the velocity components. If in addition the fluid 
incompressibility and the bottom impermeability are satisfied identically, the equations at 
the surface can be derived from Zakharov's Hamiltonian [60]. Thus, both principles naturally 
assume that the flow is exactly irrotational, as it is the case of the water wave problem formu- 
lation, but Zakharov's Hamiltonian is more constrained than Luke's Lagrangian. Luke's and 
Zakharov's variational formulations require that part or all of the equations in the bulk of the 
fluid and at the bottom are satisfied identically, while the remaining relations must be approx- 
imated [12]. It is because the irrotationality and incompressibility are mathematically easy to 
fulfill, that they are chosen to be satisfied identically. Beside simplicity, there are generally no 
reasons to fulfill irrotationality and/or incompressibility instead of the impermeability or the 
isobarity of the free surface, for example. In this paper, we illustrate the benefit of releasing 
the constraints of exact irrotationality and incompressibility, since approximations of these 
relations are sufficient in most practical cases. 

Variational formulations involving as few dependent variables as possible are often regarded 
as simpler [59]. It is understandably tempting to solve exactly (i.e., analytically) as many 
equations as possible in order to 'improve' the solution accuracy. This is not always a good 
idea, however. Indeed, numerical analysis and scientific computing know many examples 
when efficient and most used algorithms do exactly the opposite. These so-called relaxation 
methods — e.g., pseudo-compressibility for incompressible fluid flows [26] — have proven to be 
very efficient for stiff problems. When solving numerically a system of equations, the exact 
resolution of a few equations does not necessarily ensure that the overall error is reduced: 
what really matters is that the global error is minimized. A similar idea of relaxation may 
also apply to analytical approximations. 

In this study we would like to elucidate the benefit of using relaxed variational methods for 
the water wave problem. In other words, we illustrate the advantage of using a variational 
principle involving as many dependent variables as possible. We emphasize that our primary 
purpose here is to provide a generalized framework for deriving model equations for water 
waves. This methodology is explained on various examples, some of them being new to our 
knowledge. However, the potential of the present approach is far from being fully exploited. 

The paper is organized as follows. In Section 2, Luke's Lagrangian is relaxed to incorporate 
explicitly more degrees of freedom. This modification yields the Hamilton principle in its most 
general form. The advantage of this formulation is subsequently illustrated with examples 
over a fixed horizontal bottom, for the sake of simplicity. We begin in Section 3 with shallow 
water models, where some well-known and a few novel models are derived from the same 
Lagrangian, but with different subordinate constraints. In Section 4, similar model equations 
are derived in the deep water limit. In particular, a generalization of the Klein-Gordon 
equation and a remarkably simple accurate approximation for traveling waves are derived. 
A generalized ansatz, including the shallow and deep waters as limiting particular cases is 
presented in Section 5. Further generalizations are discussed in Section 6 and their advantage 
is illustrated with a variant of Serre's equations. Finally, conclusions and perspectives for 
future studies are outlined in Section 7. 

2. Generalized variational formulation 

Water wave problem possesses several variational structures [57, 36, 60]. In the present 
study we will extensively exploit the Lagrangian variational formalism. Surface gravity wave 
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equations, for a potential flow with an impermeable bottom and an impermeable free surface 
where the pressure is constant (taken to be zero), can be derived from the "stationary point" 
(point where the variation is zero) of the following functional [36]: 

£ = f 2 ! J£p d 2 xdt, J? = - T [ gy + <fr t + \{^<P? + \cfr 2 ] dy, (2.1) 
Jt x Jn J-d 

with x = (x\,X2) the horizontal Cartesian coordinates, y the upward vertical coordinate, t 
the time, V the horizontal gradient, f2 the horizontal domain, <fr the velocity potential and 
g > the acceleration due to gravity; y = r/(x,t), y = and y = —d(x,t) being, respectively, 
the equations of the free surface, of the still water level and of the bottom. For the sake of 
simplicity, the surface tension effect is neglected and the fluid density p is constant (and can 
thus be set to unity without any loss of generality), but this is not a limitation for the purpose 
of this paper. 

Integrating by parts, then neglecting the terms at the horizontal and temporal boundaries 
because they do not contribute to the functional variations (this will be done repeatedly 
below without explicit mention), Luke's variational formulation (2.1) can be rewritten with 
the following Lagrangian density: 

J2f = ht + 4>d t - \gif + \gd 2 - T [I(V0) 2 + \<fr 2 y ] dy, (2.2) 

J — d 

where the over 'tildes' and 'wedges' denote, respectively, the quantities written at the free 
surface y = rj and at the bottom y = —d. We shall also denote with 'bars' the quantities 
averaged over the water depth, e.g. 

u(x,t) = — / u(x, y, i) dy. 

( ' n(x,t) + d(x,t) J_ d{X)t) 1 1 y 

The variational formulations (2.1) and (2.2) impose that any approximation is exactly 
irrotational, i.e., the choice of an ansatz for <f> necessarily implies an irrotational motion. 
Note that the term \gd 2 in (2.2) can be omitted because, d being prescribed, it does not 
contribute to the variation process. Note also that Luke's Lagrangian (2.1) can be extended 
to the case where the bottom function d(x,t) is unknown if a condition at y = —d is added 
to the problem [54]. 

To give us more freedom while keeping an exact formulation, the variational principle is 
modified (relaxed) by introducing explicitly the horizontal velocity u = V</> and the vertical 
one v = (fry The variational formulation can thus be reformulated with the Lagrangian 
density 

JSf = fat + M " \arf - P [\{u 2 + v 2 ) + v (V0 - u) + v{<fr y - v)) dy, (2.3) 

J-d 

where the Lagrange multipliers and v have to be determined. By variations with respect 
of u and v, one finds at once the definition of the Lagrange multipliers: 

fi = u, v = v. (2-4) 

Using these definitions, (2.3) becomes 

rv 

% = <t>rj t + <j>dt - \grj 2 + / [\u 2 + \v 2 - u-Vcfr - vcfr y ]dy. (2.5) 

J-d 
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However, it is advantageous to keep the most general form of the Lagrangian (2.3), because it 
allows to choose ansatz for Lagrange multipliers /_t and v that are different from the velocity 
field u and v. Indeed, the Lagrangian density (2.3) involves six variables {rj, (f>, u, v, fx, u}, 
while the simplified Lagrangian (2.5) involves only four variables {ry, 0, u,v} and the original 
Lagrangian (2.2) only two (77 and (f>). These additional variables introduce additional freedom 
in the construction of approximations, thus allowing more subordinate relations to be fulfilled. 
The Lagrangian density (2.5) was used by Kim et al. [28] to derive the 'irrotational' Green- 
Naghdi equations for long waves in shallow water. The more general Lagrangian density (2.3) 
provides more flexibility to derive model equations, as illustrated below. 

The connection of (2.3) with the variational formulation of the classical mechanics can be 
seen applying Green's theorem to (2.3) that yields another equivalent variational formulation 
involving the Lagrangian density 

Se = (i] t + jx-Vr] -D)<f> + (dt + (i-Vd + v)4> - \grf 
rv 

+ / [fi-u-\u 2 + vv-\v 2 + (V • Li + v y )(f>]dy, (2.6) 

J-d 

which in deep water limit (d — > 00) becomes (if 4> decays faster than y _1 as y — > —00) 

Jz? = (rjt + jl • Vrj - v) (f> - i g rf 

+ / [11 • u - \u 2 + uv - \v 2 + (V • 11 + Vy) (f>] dy. 
J —00 

Thus, in the special case 11 = u and v = v , the Hamilton principle of classical mechanics is 
recovered in Eulerian description of motion, i.e., the variational principle involves the kinetic 
energy minus the potential energy plus some constraints for the fluid incompressibility, for 
the flow irrotationality and for the bottom and surface impermeabilities, as already pointed 
out by Miles [41]. In other words, the Lagrangian density (2.6) is the Hamilton principle in 
its most general form for irrotational surface gravity waves. 

Note that, via the integration by parts, the term r] t (fi, for example, can be replaced by — 77^ 
in all the Lagrangian densities given above, without loss (nor gain) of generality. Note also 
that the relaxed variational formulations involving (2.3) and (2.6) being strictly equivalent, 
one should use the more convenient one depending on the problem under consideration. 
Note finally that extensions of (2.3) and (2.6) including, e.g., obstacles, surface tensions and 
stratifications in several homogeneous layers are straightforward generalizations. For instance, 
to include the surface tension it is sufficient to add the term — cr(\/l + (V77) 2 — 1) into the 
definition of the Lagrangian density (2.6), a being the surface tension coefficient. 

The goal in this paper is to illustrate the power of the relaxed variational principle via 
some simple examples. We shall thus consider, for simplicity, a fixed horizontal bottom and 
we shall derive various approximate equations for shallow and deep waters. Possibilities for 
arbitrary depths and some other generalizations will also be discussed. 

3. Examples in shallow water 

Let consider first the shallow water case, with constant depth for the clarity of exposition. 
We introduce a realistic ansatz for these waves and then apply several constraints to derive 
various approximations, some of them being well-known, others being new. 
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3.1. Choice of a simple ansatz. For a long wave in shallow water, in potential motion on 
a horizontal impermeable sea bed at y = —d, it has long been noticed that the velocity field 
can be well approximated truncating the following expansion due to Lagrange [33]: 

u = u - \{y + d) 2 V 2 ii + ±(y + d) 4 V 4 u + ■■■ . (3.1) 

All Lagrange's followers (e.g., Airy, Boussinesq, Rayleigh and many others) used this type of 
expansions to derive their respective approximations [13]. Reviews on shallow water approx- 
imations can be found in Bona et al. [5, 6], Kirby [29], Madsen and Schaffer [39], Wu [58], 
Dougalis & Mitsotakis [16], among others. 

We consider here a simple ansatz of polynomial type, that is a zeroth-order polynomial 
in y for (j> and for u, and a first-order one for v, i.e., we approximate flows that are nearly 
uniform along the vertical direction. Our ansatz thus reads 

4> « <f}(x,t), u rj u(x,t), v w (y + d) (rj + d)~ x v(x, t). (3.2) 

Such ansatz are the basis of most shallow water approximations. We have also to introduce 
suitable ansatz for the Lagrange multiplier fi and v. Since fi = u and v = v for the exact 
solution, a natural ansatz for the multipliers is 

fi Ri p,(x,i), v rj (y + d) (rj + cTp 1 u(x, t). (3.3) 

With the ansatz (3.2) and (3.3), the Lagrangian density (2.6) becomes 

S£ = {r]t + iJ,-Vr,)4> - \gr, 2 

+ (ri + d)[p,-u - \u 2 + \Dv - \v 2 + 4>V -fi] . (3.4) 

Using the Green formula, the variational problem can also be written such that the Lagrangian 
density is in the following simpler form 

££ = 4>rit - \gr] 2 + (77 + d) [ p, • u - \u 2 + \ vv - \v 2 - fx- V$] . (3.5) 

The two Lagrangian densities (3.4) and (3.5) differing by a divergence term, they yield exactly 
the same equations. Thus, depending on the constraints, we use the Lagrangian density 
leading to the simpler expression. We now investigate equations led by this shallow water 
model under various subordinate relations. 

3.2. Unconstrained approximation. Without further constraints, the Euler-Lagrange of 
(3.5) yield 

(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 

\v 2 - fi- V0 - t - gr,. (3.11) 

The relations (3.6)-(3.9) imply that the motion is exactly irrotational, but the fluid incom- 
pressibility is not satisfied identically. With these four relations, the last two equations can 
be rewritten in the form: 

ht + V- [hu] = 0, (3.12) 
u t + \ Vu 2 + gVh = 0, (3.13) 



5 u : 
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- u, 


5v : 
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- v, 


5 ji : 
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5D : 
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54> : 
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5 rj : 
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where h = rj + d is the total water depth. Equations (3.12)-(3.13) are the very well-known 
nonlinear shallow water equations, also known as Airy or Saint- Venant equations (Wehausen 
& Laitone [56], §28). They are sometimes called non-dispersive fully-nonlinear approxima- 
tion because their classical derivation assumes long waves without the extra hypothesis of 
small amplitudes. These equations have a canonical Hamiltonian structure for the conjugate 
variables rj and <p with the Hamiltonian 

f^{\gr? + Hv + d) (V^) 2 } d 2 x. 

The Saint- Venant equations do not admit smooth progressive wave solutions. They are 
nevertheless widely used because they can be solved analytically by the method of charac- 
teristics [50]. Moreover, numerous efficient finite volumes type schemes have been proposed 
[62]. These equations have also been validated by comparisons with experimental data even 
for the wave runup case [52]. Below we present several more elaborated models. 

3.3. Constraining with free surface impermeability. We now constrain the ansatz (3.2) 
imposing that the impermeability of the free surface is satisfied identically. Since the surface 
impermeability is expressed through the velocity (/J-,v) in (2.6), we substitute 



v = rj t + A-Vr/, (3.14) 

into the Lagrangian density (3.4), and the subsequent Euler-Lagrange equations are 

5u : = p — u, (3.15) 

5v : = rjt + p-Vr] — v, (3.16) 

5p: = u + \vVn - V0, (3.17) 

54>: = rjt + V ■ [{ri + d) fi} , (3.18) 

5n: = p-u — \v? — \v 2 — p • Vi^ — 4>t — 9V 

- l(r] + d)[v t + p-Vv + vV-p}. (3.19) 



The relations (3.15) and (3.17) link the velocity potential and the horizontal velocity as 
7^ u = p and, therefore, equations (3.15)-(3.19) cannot be derived from Luke's variational 
principle. Relations (3.15) and (3.18) provide the mass conservation and hence, with (3.14), 
the approximation (3.15)-(3.19) implies that the fluid incompressibility is fulfilled identically. 

Eliminating (ft, p and v from the horizontal gradient of (3.19), the system (3.15)-(3.19) 
becomes 

ht + V • [hit] = 0, (3.20) 
u t + \Vu 2 + gVh + \h~ l V[/i 2 7] = (u • Vh) V(W • u) - [u • V(/iV • u)} Vh, (3.21) 

with h = n + d and where 

7 = v t + u ■ Vv = h { (V • u) 2 - V • u t - u • V [ V • u] } , (3.22) 

is the fluid vertical acceleration at the free surface. 

In the two-dimensional case (one horizontal dimension) the right-hand side of (3.21) van- 
ishes and the system (3.20), (3.21) reduces to the equations first derived by Serre [49], inde- 
pendently rediscovered by Su and Gardner [51] and again by Green, Laws and Naghdi [22]. 
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It is sometimes called weakly- dispersive fully-nonlinear approximation [58]. These equations 
admit a traveling solitary wave solution 

77 = a sech 2 \x(x\ — ct), c 2 = g(d + a), (xd) 2 = 3a(d + a) _1 , 

which is linearly stable [35]. Note that this solution does not impose any limitation on the 
wave amplitude, meaning that Serre's equations are physically inconsistent for the highest 
waves. Note also that the Serre equations have a non-canonical Hamiltonian structure [35]. 

In three dimensions, equations (3.20)-(3.21) were called by Kim et al. [28] 'irrotational' 
Green-Naghdi equations. If the right-hand side of (3.21) is neglected, we recover the classical 
Green-Naghdi equations [23]. 

Craig and Grooves [11], and many others, also derived various shallow water models starting 
from a variational principle. The main difference with our approach is that we do not introduce 
any small parameter but we make assumptions on the vertical structure of the flow. 

3.4. Constraining with incompressibility and partial potential flow I. Here, we re- 
strict the freedom imposing that the velocity potential is related to the horizontal velocity as 
u = V<^>, and that the fluid incompressibility V • u + v y = is fulfilled, together with the 
relations /x = u and v = v, i.e., we take the subordinate conditions 

ft = u, v = v, u = V</>, v = — (77 + <i) V 2 0. 

These constraints do not impose exact irrotationality because v 7^ <p y . Obviously, we shall 
derive an approximation which lies "between" the Saint- Venant and Serre equations. 
Thus, the Lagrangian density (3.5) becomes 

JS? = 4>m - \gr, 2 - i(r / + d)(V0) 2 + \^ + df{V 2 4>) 2 , (3.23) 

and its Euler-Lagrange equations yield 

50: = 7] t + V- [(v + d) V0] + fV 2 [(j) + rf) 3 (V 2 <^)] , 

S V : = & + g V + i (V0) 2 - + d) 2 (v 2 4>) 2 . 

It seems that these equations have never appeared before in the literature. They are a 
generalization of the so-called Kaup-Boussinesq (or canonical Boussinesq) equations [27, 32] 
and are thus referred to as the gKB equations. This can be seen noticing that the gKB 
equations can be derived from the canonical Hamiltonian 

J { \ g V 2 + \ ( V + d) (V^) 2 - i (r, + d) 3 (V 2 0) 2 } d 2 z, (3.24) 

while the classical Kaup-Boussinesq (cKB) equations are obtained replacing (77 + d) 3 by d 3 in 
(3.24) and restricting the resulting Hamiltonian to one horizontal dimension. Note that the 
Lagrangian (3.23) is cubic in 77, meaning that it has (at most) one local minimum in rj and 
not a global one. 

The linearized gKB and cKB systems admit the special traveling wave solution 

77 = acos k{x\ — ct), c 2 = gd(l — ^k 2 d 2 ), (3.25) 

implying that these equations are linearly ill-conditioned (c 2 < for kd > y/3). However, if, 
like the cKB, the gKB equations are integrable, they may be a somewhat interesting model 
for gravity waves in shallow water. 
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3.5. Constraining with incompressibility and partial potential flow II. So far, all 
the approximations derived turned out to be such that n = u and v = v. We propose here a 
novel approximation that does not satisfies one of these identities and that is an interesting 
variant of the previous model. 

We impose a partially potential flow such that fi = V</> and v = cp y , together with the 
incompressibility condition V • u + v y = and the condition u = /x. Thus, substituting the 
constraints 

fi = u = V</>, v = 0, v = -(r) + d)V 2 4>, 
into the Lagrangian density (3.5) yields 

2 = 4>rj t - \gr, 2 - l( v + d){V^) 2 - ± ( V + df (V 2 ^) 2 , 

and the corresponding Euler-Lagrange equations are 

54> : = 7] t + V- [(v + d) V0] - | V 2 [(ri + d) 3 (V 2 0) ] , (3.26) 

S V : = fa + gri + | (V0) 2 + i(r? + d) 2 (V 2 0) 2 . (3.27) 

These equations can be derived from the canonical Hamiltonian 

J Q {hav 2 + \{v + d)(V4>f + \{n + df (v 2 ^) 2 } d 2 *, 

which is always positive (an interesting feature for modeling water waves). To the linear 
approximation, equations (3.26), (3.27) have the progressive wave solution 

j] = acosk(xi-ct), c 2 = gd(l + \k 2 d 2 ), (3.28) 

which is well-behaved (i.e., c 2 is never negative). Comparisons with the gKB equations suggest 
to refer to equations (3.26)-(3.27) as regularized general Kaup-Boussinesq (rgKB). However, 
the linear dispersion relation (3.28) approximates the dispersion relation of linear waves - 
i.e., c 2 = g tanh(kd) / k — only to the order 0(k 2 ), while (3.25) is 0(A: 4 ). Therefore, the rgKB 
equations are not very interesting for modeling water waves, but these equations may be of 
interest to model other physical processes. 

3.6. Other constraints and generalizations. We briefly discuss now some other possible 
generalizations . 

3.6.1. Constraining with incompressibility and potential flow I. In the previous example, we 
have constructed an approximation such that [i = u but v ^ v. Now, we release the constraint 
fi = u and keep the other constraints. Thus, we impose 

pi = Vcj), v = 0, v = — (rj + d) V • u, 

so that the pseudo velocity field (/x, v) is irrotational while the velocity field (u, v) is incom- 
pressible. After same elementary algebra, the Lagrangian density becomes 

J2f = 4>h t - \gn 2 + hu-Vj> - \hu 2 - \ h 3 (V • uf - /i(V0) 2 , 

where h = rj + d. The Euler-Lagrange equations give 

Su: = hV4> - hu + ±V [/i 3 V • u] , 

54>: = h t - V-[hu] + 2V-[/tV0], 

Srj: = 4> t + gr] + \u 2 + (V^) 2 - u ■ + \ h 2 (V • u) 2 . 
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The linearization of this system of equations have a (27r//c)-periodic sinusoidal traveling wave 
solution with the dispersion relation 

c 2 = ^(l + ffcVKl + ffcV)- 1 = gd(l + \k 2 d 2 ) + 0(k 4 ), 

which, like the previous example, is not satisfactory for water waves. However, these equations 
may be of interest in other contexts than water waves. 

3.6.2. Constraining with incompressibility and potential flow II. We now assume that the 
pseudo velocity field (/x, v) is divergence free, while the velocity field (u, v) is irrotational, 
i.e., we impose the constraints 

u = V0, v = 0, v = — (77 + d) V • fi. 

The Lagrangian density becomes 

jSf = 4> m - \gtf - i( v + d){V4>) 2 , 

which yields the Saint- Venant equations. Thus, these constraints do not bring anything new. 
It should be emphasized that this is the case for the special shallow water ansatz we are 
considering here, but this is not necessarily the case for other ansatz. 

3.6.3. Further possibilities. The constraints of Sections (3.4) to (3.6.2) can be unified into 
a single formalism considering combinations. Indeed, the velocity field (u, v) being not 
more (nor less) physical than the pseudo-velocity field (n, v) and the potential velocity field 
(^4>,4> y ), the constraints can be imposed by combinations of these three fields. For instance, 
we could impose the irrotationality for the field 

( C\U + C 2 jU + (1 — Cl — C 2 ) V(f) , C\V + C 2 f + (1 — Cl — C2)(j)y ), 

the fluid incompressibility for the field 

( C 3 U + C4/X + (1 - C 3 - C 4 ) , C 3 1> + C^U + (1 - C 3 - C4)(f)y ), 

and so on for any constraint we may think of. The c n are parameters at our disposal. We can 
choose them in a convenient way based on some mathematical and physical considerations. 
For example, imposing that the approximate equations derived must be linearly well-posed 
and/or have better dispersion relation properties. 

In the examples above, only some kinematic constraints (irrotationality, incompressibility, 
impermeability) were used. We could have also considered dynamical constraints based on, 
e.g., the Bernoulli equation, or other relevant dynamical equations. 

The relaxed variational principle provided a common platform for deriving several shallow 
water equations from the same ansatz in changing only the constraints. Beside the ansatz, no 
further approximation were made and the derivations required only some elementary algebra. 
Using more general ansatz — i.e., involving more free functions and parameters — one can 
introduce more constraints, if desired, and derives more accurate approximations. A simple 
example is given in Section 6 below. 

4. Examples in deep water 

We illustrate here the advantages of the relaxed variational principle in the opposite limiting 
case of deep water. 
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4.1. Choice of an ansatz. For progressive waves in deep water, the Stokes expansion shows 
that the velocity field varies nearly exponentially along the vertical (Appendix A). Even for 
very large unsteady waves (including breaking waves), accurate numerical simulations and 
experiments have shown that the vertical variation of the velocity field is indeed very close to 
an exponential [24, 25]. Thus, this property is exploited here to derive simple approximations 
for waves in deep water. 

Let k > be a characteristic wavenumber corresponding, for example, to the carrier wave 
of a modulated wave group or to the peak frequency of a JONSWAP spectrum. Following 
the discussion above, it is natural to seek approximations in the form 

{0; u;v; fi; v} « {fi; u; v ; ft; v} e re ^, (4.1) 

where (ft, u, v, ft and v are functions of x and t that will be determined using the variational 
principle. The ansatz (4.1) is certainly the simplest possible that is consistent with experi- 
mental evidences. This ansatz has already been used by Kraenkel et al. [31] for building their 
approximation. Possible generalizations are discussed in Section 6. 

For the sake of simplicity, we introduce the constraints ft = u and v = v. Thus, the ansatz 
(4.1) substituted into the Lagrangian density (2.5) yields 

2fti? = 2Kcpr/ t - gnrj 2 + \u 2 + \v 2 - u • ( - KcjyVr]) - KV(f>. (4.2) 

With (or without) subordinate relations, this Lagrangian gives various equations. We inves- 
tigate two cases here. 

4.2. Unconstrained approximation. Without further constraints, the Euler-Lagrange equa- 
tions yield 

5u : = u — + Kcp'Vr], 
5v: = v — K(p, 

5 4> ■ = 2Kr]t + 'V • ii — kv + ku • Vr/, 

5 7] : = 2gn rj + 2k fit + k V • ( ^ u ) . 

The two first relations imply that this approximation is exactly irrotational and their use in 
the last two equations gives 



:4> [V 2 rj + k(Vt]) 2 ' 
- k4> 2 Vt] 



(4.3) 
(4.4) 



Since equations (4.3)— (4.4) derive from an irrotational motion, they can also be obtained 
from Luke's Lagrangian (2.2) under ansatz (4.1). Equations (4.3)— (4.4) are a deep water 
counterpart of Saint- Venant equations for shallow water waves; this claim will appear clearer 
in Section 5. They can also be derived from the canonical Hamiltonian 



hv 2 + 1« 1 



•Vrj 



d 2 x. 



(4.5) 



This 'simple' Hamiltonian is quartic in nonlinearities and involves only first-order derivatives. 
It has to be compared with Zakharov's quartic Hamiltonian (B.2) which involves second-order 
derivatives and pseudo-differential operators. However, Zakharov's quartic Hamiltonian is 
valid for broad spectra. Note that the Hamiltonian (4.5) cannot be derived from the exact 
one (B.l), since the latter assumes that irrotationality and incompressibility are both satisfied 
identically in the bulk, while the incompressibility is not fulfilled by equations (4.3)-(4.4). 
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To the linear approximation, after elimination of eft, equations (4.3)-(4.4) yield 

Vtt ~ i(5/«)V 2 r? + \gnn = 0, (4.6) 

that is a Klein-Gordon equation. For this reason, equations (4.3)-(4.4) will be referred here as 
generalized Klein-Gordon (gKG). The Klein-Gordon equation is prominent in mathematical 
physics and appears, e.g., as a relativistic generalization of the Schrodinger equation. The 
Klein-Gordon equation (4.6) admits a special (2-7r/A;)-periodic traveling wave solution 

n = a cos k(xi — ct), c 2 = \ g (k 2 + k 2 ) (nk 2 )^ 1 . 

Therefore, if k = k the exact dispersion relation of linear waves (i.e., c 2 = g/k) is recovered, 
as it should be. This means, in particular, that the gKG model is valid for spectra narrow- 
banded around the wavenumber k. 

We focus now on (2-7r/K)-periodic progressive waves solution of the gKG equations, i.e., we 
seek for solutions depending only on the variable 6 = K,(x\ — cb). We were not able to find 
an exact analytic solution but a Stokes- like expansion gives some interesting insights. To the 
seventh-order, we have 

Kr] = acos6 + ia 2 (l + ffa 2 + ^^a 4 ) cos26> 

+ | a 3( 1 + 99 a 2 + 11807 a 4) cos 3^ + ^ a ^ 1 + f a ^ cosAd 

+ gfa 5 (l + ^a 2 ) cos 56 + §a 6 cos6# + ga 7 cos 76 + 0(a 8 ), 

g-ljj, = a (l_l a 2_|9 a 4_^ a 6^ sin ^ + l a 2^ + ll a 2 + 547 a 4^ sin20 

+ l« 3 ( 1 + ^« 2 + W« 4 ) sin30 + ^(i + lfa 2 ) sm 40 

+ ifa 5 (l + ^fa 2 ) sin 56 + fa 6 sin 66 + ga 7 sin 70 + 0(a 8 ), 

g -\ K \c = 1+ ia 2 + i a 4 + ffa 6 + 0(a 8 ). 

The expansions of n and 4> match the exact Stokes wave (c.f. Appendix A) up to the third- 
order (non-matching coefficients are displayed bold). This is not surprising since the gKG 
equations are cubic in nonlinearities. A bit more surprising is that the phase velocity c is 
correct up to the fifth-order. But the most interesting is that, to the leading order, the n-th 
Fourier coefficient is (for all n up to infinity) 



2™- 1 (n-l)!' v ' ; 

which is also the case for the exact Stokes wave (Appendix A). 

In comparison, for the cubic Zakharov equations (B.3)-(B.4), the phase velocity is correct 
only up to the third-order and the Fourier coefficients do not verify the asymptotic behavior 
(4.7) (see Appendix B). Truncating Zakharov's Hamiltonian at the order N + 1 in nonlin- 
earities, the corresponding Stokes double series is correct up to the order N in the expansion 
parameter. But none of these higher approximations have the exact asymptotic behavior 
(4.7) for their Fourier coefficients because they involve expansions around n = 0, while the 
gKG does not. Zakharov's equations are valid for broad spectra, however, unlike the gKG 
equations. 

4.3. Constraining with the free surface impermeability. In order to satisfy the free 
surface impermeability identically, we take 

v = rj t + it • V77, 



5u : 





5<f> : 





5 r] : 
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and the Lagrangian density (4.2) becomes 

2k^£ = 4>(KT] t + V • u) - gKT] 2 + \u 2 + l(ryt + u • Vrj) 2 , (4.8) 

while the Euler-Lagrange equations yield the relations 

u + (r] t + u- Vrj) Vr, - V0, (4.9) 

Krj t + V • u, (4.10) 

2g K rj + K<j) t + rj tt + (u ■ Vrj) t + V • (urj t ) + V • [{u ■ Vrj)u]. (4.11) 

The relation (4.10) implying that V • u + v y = 0, the solution satisfies the incompressibility 
identically. On the other hand, the irrotationality being not verified identically, equations 
(4.9)-(4.11) cannot be derived from Luke's variational formulation. Note that (4.9) yields 
V(/> = u + v'Vr] that is exact for potential flows [12, 20]. 

As for the shallow water case, the potential 4> can be eliminated from equations (4.9), thus 
yielding a deep water analog of Serre's and the Green-Naghdi equations. 

To the linear approximation, relations (4.9)-(4.11) can be combined into a single equation 
for the elevation of the free surface: 

(V 2 -K 2 ) Vtt + 2gnV 2 r, = 0, 

which admits the special (2-7r//c)-periodic solution 

r\ = a cos k(x\ — ct), c 2 = 2 g k (k 2 + k 2 )^ 1 . 

Therefore, if k = k the exact linear approximation is recovered, as it should be. Again, this 
means that this model is valid for narrow-banded spectra. How narrow will be investigated 
now. 

4.3.1. Two-dimensional progressive waves. We seek now exact solutions depending only on 
£ = xi — ct (two-dimensional progressive waves). Equations (4.9)-(4.10) yield respectively 

4>£ = u + (u — c) r/£ , u = cur] + cK\, 

where K\ is an integration constant. Substituting these relations into (4.11) and multiplying 
the result by n^, after one integration and some algebra, we obtain 

'drA 2 K 2 - (nri - K\ s~ 2 ) 2 



dU e- 2 {Ki] + K x - l) 2 



where 



e 2 = 2g K - 1 c~ 2 - 1, 



(4.12) 



and where Kq is another integration constant. Assuming that ??g = at the wave crest where 
rj = a and at the trough where r] = —b (hence a + b is the total wave height), we find 

Kq = ^K(a + b), K\ = ^K(a — b)e 2 . 

Solving equation (4.12), together with the condition of zero mean free surface elevation, we 
obtain an exact (27r/fc)-periodic solution in the parametric form [21]: 

&£(t) = t — ee~ l sinr, k t]{t) = \ e 2 e^ 1 + e cost, (4.13) 

with 

e = \k(a + b), fctfT 1 = e + \ (e _1 +e)e 2 , 
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e being a wave steepness and r being a parametric variable such that a crest is at r = 
and the closest troughs at r = ±7r. This remarkably simple solution describes surface waves 
as trochoids. Smooth surface profiles are obtained for ^ e < e, e <C e corresponding to 
quasi-sinusoidal solutions (infinitesimal waves). For the limiting case e = e, in the vicinity of 
the crest k(a—rj) ~ ^e(6/c£)3 so the solution involves a sharp angle forming a 0° inner angle 
(i.e., a cusp), while the exact angle should be 120° for irrotational traveling waves. 
So far, k is a free parameter at our disposal. We shall now investigate various choices. 

4.3.2. Simple approximation. A 'natural' choice is to take k = k, yielding the steepness 
e ~ 0.596 for the limiting wave. With this peculiar choice of k, a Stokes-like expansion of the 
solution of (4.13) is 



k rj 



1^2, 



acosfe; + T^a 



l + lfa 2 + f§fa 4 ) cos 2^ 



+ |a 3 l 



l + ffa 2 + ^a 4 ) cos3fcf + |a 4 (l + f§ 



+ j||a 5 (l + f|a 2 ) cos5/c£ + §£a 6 cos6/c£ + 



a ' ) cos 4/c£ 
§g£a 7 cos7A£ + 0(a 8 ), 



1 + \o? + ia 4 



+ §±a 6 + 0(a s 



As for the gKB approximation, the Fourier coefficients of this Stokes-like expansion satisfy the 
asymptotic expression (4.7) and the expansion of c is exact up to the fifth-order. Therefore, 
taking k = k leads to a quite accurate approximation. However, another value of k may lead 
to an even more interesting approximation. 

4.3.3. Improved parameter. With the trochoidal solution and the relation between the pa- 
rameters obtained above, the Lagrangian density (4.8) integrated over one wavelength yields 
the Lagrangian 

2 



4fc 3 



TTj? 
k 3 



1 + 



K C 



25 



+ 



k 2 , 



2 9 k 4 5 2 



(4.14) 



Thus, the Lagrangian is minimum when the steepness is minimum or, equivalently, the wave- 
length is maximum if the wave height is kept constant. The variational principle being defined 
with fixed horizontal and temporal boundaries, an improved parameter n can be obtained min- 
imizing L keeping k and c constant. Thus, after some algebra, the equation dL/dn = gives 
two possible solutions for k: 



k± = V2k 1 + 2e 2 ± VI - 4e 2 



(4.15) 



provided that e ^ 1/2. Both solutions correspond to an extremum of £. However, a solution 
is stable only if k is a minimum of L, i.e., if 

d 2 £ 



dK 2 



> o, 



K,=K,± 



which, after some algebra, yields the condition 

v^ £ K± < k. 

The first solution 



(4.16) 



k + I k = V2 1 + 2e 2 + \A - 4e 2 



1 + i e 4 + ( £ 6) 
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fi(0)/c 



FIGURE 1. Deep water phase velocity as a function of the fluid speed at crest for the 
improved solution derived in §^.3.3. Solid line: n+ solution; Dashed line: k- solution. 



is very close to the wavenumber k when the steepness e is small. When the steepness in- 
creases form to 1/2, the dimensionless parameters n/k, kc 2 /g, k 3 L/g and \/3£K/k in- 
crease monotonically. The highest wave is obtained for e = 1/2, where n/k = 2/\/3 ~ 1.15, 
kc 2 /g = 3\/3/4 « 1.3 and V3eK/k = 1. All these waves are smooth and should be stable 
because the condition (4.16) is fulfilled. 
The second solution 



/ k = V2 \ 1 + 2e 2 - y/l - 4e 2 



is admissible (i.e., real) if 1/2 ^ e ^ \J \/5 — 2 ~ 0.486. As the steepness decays from 1/2 to 
Vv5 — 2, the crest sharpens, the limiting value e = e = y \/5 — 2 corresponding to a cusp 
at the crest. Thus the sharp-crested wave is not the highest one. All these waves should be 
unstable because the condition (4.16) is violated, i.e., the solution k_ is a maximum of the 
Lagrangian and not a minimum. 

Note that the relations linking the parameters can be conveniently rewritten 



kc 2 2e £2 ^ ^ {k[n±) + ^/(k/ K± ) 2 -2e 2 -i* 

~~i7F + ~' £ ~ 2~T7 2 1 

which, together with (4.15), allow to express kc 2 /g as function of e only. Another parameter 
of physical interest is the fluid velocity at crest £t(0) divided by the phase velocity c, i.e., 

u(0) _ e (e + 2e + ee 2 ) k± 
~c~ ~ 2ik ' 

The variations of c 2 with respect to this parameter are displayed on Figure 1. 

4.3.4. Remarks. Taking the improved parameter k has increased somewhat the accuracy of 
the approximation; the main gain is qualitative, however. Indeed, we found that waves 
should be unstable before a sharp crest is formed. Such behavior is not predicted by low- 
order perturbation expansions. This simple example is a remarkable illustration of the power 
of the variational method. The trochoidal wave described here is probably the approximation 
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with the highest ratio accuracy / complexity ever derived for a traveling wave in deep water. 
This approximation has first been derived by Geniet [21] via a different approach. 

5. Arbitrary depth 

A general ansatz, for waves in finite constant depth and satisfying identically the bottom 
impermeability, is suggested by the linear theory of water waves: 



cosh kY 
cosh nh 



(x,t), 



u 



cosh kY 

cosh nh 
cosh kY 



-u(x, t) 



sinh kY 

sinh nh 
sinh kY 

sinh nh 



v(x,t), 



(5.1) 



cosh nh 1 

where Y = y + d and h = rj + d. The parameter n is a characteristic wave number to be 
made precise a posteriori. This ansatz is uniformly valid for all depths because it yields the 
shallow water one (3.2) as k — > 0, and the deep water one (4.1) as d — > oo. Obviously, the 
ansatz (5.1) is valid for wave fields with wavenumber spectra that are narrow-banded around 
K. Substituting the ansatz (5.1) into the relaxed variational principle (2.6), we obtain 

sinh(2/c/t) — 2nh 



= [rjt + ft- Vtj] 



hn 2 



+ 



v v 



+ [A 



u 



+ b 



\u 2 + 



2nh 



)V • [i — Ktanh(K/i)0 jl • Vt? 



2KCOsh(2K/i) — 2k 

sinh(2K/i) + 2nh 



2KCosh(2K/i) + 2k 



sinh(2K/t) ' ^ 5 ' 2 ^ 

Applying various constraints, one obtains generalized equations including the ones derived 
in Sections 3 and 4 as limiting cases. In particular, we can derive arbitrary depth general- 
izations of the Saint- Venant and 'irrotational' Green-Naghdi equations. Thus, the relaxed 
variational principle leads simple derivations of model equations for finite amplitude waves in 
arbitrary depth, which is not the case with the classical perturbations techniques. These de- 
velopments are left to future investigations. We note also some similarity between the ansatz 
(5.1) and the coupled-mode theory for small-amplitude water waves [2]. 

6. Generalizations 

The ansatz (5.1) can be generalized in many relevant ways, depending on the problem under 
consideration. Natural generalizations can be based on the velocity vertical variations given 
by, e.g., higher-order deep and shallow water theories, or obtained fitting some experimental 
data. In this section we propose a possible generalization based on an ansatz of the form 



cosh kY 
cosh kK 



x,t) 



u 



cosh kY 
cosh kK 



u(x,t), 



sinh kY 
sinh nh 



v(x,t), (6.1) 



where A is a parameter at our disposal. If A = 1, the ansatz (5.1) is recovered, but the 
case A / 1 does not correspond to the vertical profile predicted by any theory based on 
perturbation expansions. Still, this type of ansatz is of some interest, as we shall see below. 
Note first that in the deep water limit d — > oo, the ansatz (6.1) becomes 



e XK( - y -^4>(x,t), u » e XK{lJ -' n) u(x,t), v « e x ^ y ~ v) v(x,t), 



and thus, via the change of parameter \k i— > k, the ansatz (4.1) is recovered. This means 
that (6.1) is not more general than (5.1) in deep water. On the contrary, these two ansatz 



16 



D. CLAMOND AND D. DUTYKH 



are very different in finite depth. We illustrate this claim in the simple case of shallow water 
(nd — > 0) when (6.1) becomes 

A 



(x,t), U Ri u(x,t), V 



y + d 

rj + d 



v(x,t), 



where we have replaced 4> an d u by and u, respectively, since they are equal in this limiting 
case. 

6.1. Modified Serre's equations. For the sake of simplicity, we consider here only one 
horizontal dimension, say x±, and we set x\ = x and u\ = u, for brevity. We also consider 
the special case \x = u and v = v together with the constraint v = rjt + urj x (free surface's 
impermeability). Thus, the Lagrangian density (2.6) yields 

& = (Vt + [(V + d)u] x ) 4> - \gr} 2 + i(7? + d)u 2 + i/3 2 (r? + d) [r/ t + ut] x } 2 , (6.2) 

where f3 2 = (2A + 1) . After some algebra, the Euler-Lagrange equations leads to the 
following equations 

h t + [hu] x = 0, (6.3) 

u t + uu x + gh x + fhT x [h 2 ^] x = 0, (6.4) 

where 7 is defined in (3.22). If /? = l/\/3 the classical Serre equations are recovered. 
Equations (6.3)-(6.4) admit a (27r//c)-periodic cnoidal traveling wave solution: 

with ^ = c — ct, dn and sn being elliptic functions of Jacobi of parameter m (0 ^ m ^ 1), 
and where X = K(m) and = E(m) are the complete ellipic integrals of the first and second 
kinds, respectively (Abramowitz & Stegun 1965, #17.3, [1]). The parameter k is a sort of 
wavenumber, a is the wave amplitude (mean level to crest elevation), H is the total wave 
height (trough to crest elevation) and c is the wave phase velocity observed in the frame of 
reference without mean flow. The wave parameters are related via the relations 

vrx maK 2 gH 

k = Wr?' H = 1? p' ( xd > = a2 2 ' ( 6 - 6 ) 

IK K — h, mp A cr 



it 



gH(d + a) (d + a-H) 
g (d + a) 2 (d + a — H) — d 2 c 2 ' 1 j 

In the limiting case m — >■ 1, we have X — >• 00, — )■ 0, k — > 0, H — > a and hence, the 
classical solitary wave solution is recovered 

1 2 1 / \ 2 i 1 \ a (fixd) 2 

i] = asech ^x{x — ct), c = g{d + a), — 



d 1 - {f3xd) 2 ' 

At this stage, j3 is still a free parameter. An suitable expression for this parameter may 
be obtained substituting the solution (6.5) into the Lagrangian density (6.2), integrating 
J2? over one wavelength, then solving d£/d/3 = keeping k and c constant (as well as g 
and d), the other parameters varying according to relations (6.6-6.7). Thus, after some 
cumbersome algebra, we found f3 = for this parameter, which is not very interesting for 
practical applications. A possible alternative here is to choose f3 such that the exact relation 
c 2 = gtan(xd)/x is satisfied identically or up to some asymptotic order. 
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7. Discussion 

In this study, we have illustrated the advantage of using a variational principle with as 
many variables as possible. We call it the relaxed variational principle, since the Lagrangian 
density (2.6) involves more degrees of freedom (i.e., the variables rj, (ft, u, v, and v) com- 
pared to the two degrees of freedom (77 and (ft) in the classical case. In particular, these 
extra variables can be used to impose various constraints such as incompressibility, irrota- 
tionality, impermeability, etc. The practical use of the relaxed formulation was illustrated 
on numerous examples in shallow, deep and intermediate waters. Thus, we obtained several 
approximations, some well-known, some new to our knowledge. 

In the shallow water regime, we have first obtained the classical nonlinear shallow water (or 
Saint- Venant) equations (3.12)-(3.13). Then, with the same ansatz (3.2) but imposing the 
constraint of the free surface impermeability, we have derived the irrotational Green-Naghdi 
equations (3.20)-(3.21). Applying the incompressibility constraint and choosing differently 
the pseudo- velocity field, we have obtained two kinds of generalized Kaup-Boussinesq equa- 
tions. Several ways of further generalizations were also outlined. 

In deep water, two models were considered. Namely, we derived deep water counterparts 
of the celebrated Saint- Venant and Serre equations. The former has a canonical Hamilton- 
ian formulation and degenerates to the Klein-Gordon equation in the linear approximation; 
we thus called the new system (4.3)-(4.4) generalized Klein-Gordon equations. The latter 
could be solved analytically for a two-dimensional traveling wave. This solution is a striking 
illustration of the power of the variational formulation compared to asymptotic expansion 
methods, especially for large amplitudes when the expansion parameter is no longer small. In 
addition, both equations were shown to possess excellent asymptotic properties with respect 
to Stokes-like expansions. 

The case of arbitrary depth has also been briefly considered. In particular, it has been 
shown how easily one can introduce an ansatz valid for all depths. Indeed, the vertical 
variation of the velocity field suggested by the linear theory provides at once such a general 
ansatz which degenerates to previous cases when the water is shallow (nd —> 0) or deep 
(d — > 00). This simplicity and flexibility of the variational principle is quite remarkable 
compared to perturbation methods. 

Finally, we considered a generalized ansatz and we illustrated its consequence in the limiting 
case of shallow water. In this way, we derived the modified Serre equations and subsequently 
obtained exact cnoidal and solitary wave solutions. The main purpose of this example was to 
illustrate the fact that one can introduce an ansatz which is not inspired by any asymptotic 
expansion and nevertheless lead to reasonable approximations. 

In the present paper, some further possibilities for generalizations are also mentioned. 
However, we have to emphasize that not all ansatz and constraints will necessarily lead 
to physically relevant and tractable approximations; the same is true for models derived 
from asymptotic expansions, however. Nonetheless, the relaxed variational formulation is 
sufficiently versatile to allow easy derivations of physically sound models. We have illustrated 
this claim, in particular, by showing how it is simple to obtain approximate equations valid 
for all depths. 

Sometimes, the choice of the constraints may seem to be rather ad hoc, but that should not 
be surprising. Indeed, the water wave theory already knows several ad hoc 'tricks' intended 
to improve the approximation quality. For instance, it was proposed in [38] to replace the 
polynomial shallow water expansion (3.1) by a (m,n)-Pade approximation, the orders m and 
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n being chosen to improve the linear dispersion relation of progressive waves. Another exam- 
ple is the use of the velocity potential defined at some depth yo [43] and, as before, the free 
parameter yo is chosen to improve linear dispersion characteristics. The approach proposed 
here is not more ad hoc that any example mentioned above. Moreover, the variational princi- 
ple allows for greater flexibility in the choice of ansatz. Thus, the approximations derived via 
the relaxed variational procedure must be studied a posteriori in order to verify their math- 
ematical consistence and their physical relevance. This is also the case for approximations 
derived via perturbation techniques, and many such approximations commonly used have not 
yet been justified on a rigorous mathematical basis. 

In the several examples presented here, the ansatz involve free parameters that we have 
chosen constant for simplicity. One may also consider these parameters as functions and 
find their values from the stationary point of the Lagragian. Doing so will lead to more 
complicated equations, but this is not a major issue if these equations are intended to be 
treated numerically. 

In order to derive approximate models, variational formulations are attractive alternative 
to asymptotic expansions. However, both approaches can be also combined. Indeed, once 
the variational principle has been applied to an ansatz, asymptotic expansions can be further 
applied to obtain simpler models. For instance, one could consider 'unidirectionalized' ap- 
proximations [44, 45] to derive variants of Korteweg and de Vries [30], Dysthe [17], Camassa 
and Holm [9], Degasperis and Procesi [14], Kraenkel et al. [31], and other equations. This 
possibility will be investigated in future works. 

For the sake of simplicity, we have considered only gravity waves propagating at the surface 
of a single layer of a homogeneous fluid with a horizontal bottom. It is trivial to introduce 
a relaxed variational formulation including, e.g., surface tension, stratifications in several ho- 
mogeneous layers and obstacles. Such general variational formulations, together with relevant 
ansatz and well chosen constraints, will easily lead to interesting models. For perfect fluids, 
variational formulations can also be obtained for rotational motions [18, 36, 34, 48, 42, 10]. 
A relaxed version of such variational principles will facilitate the derivation of approximate 
models. 

The numerical models for simulating water waves are undergoing constant improvements. 
The state of the art can be found in recent reviews (e.g., [19, 15, 37]). Certainly, the varia- 
tional principle is the tool of choice to derive efficient approximations. By efficient, we mean 
models that capture most of the relevant physics and which, in the same time, can be easily 
and rapidly solved numerically. One interesting direction for future researches is the devel- 
opment of numerical schemes preserving hamiltonian structure at the discrete level. Another 
interesting application is the derivation of new improved models with uneven bathymetry for 
coastal hydrodynamics and tsunami wave modeling (Synolakis & Bernard 2006, [53]). 

The use of a variational principle for modeling surface waves is by no mean new. However, 
its power has not yet been fully exploited. The present paper is a further contribution in this 
direction. Obviously, advantages of relaxed variational principles may as well be exploited in 
physical contexts other than water waves, such as plasma physics, nonlinear optics, etc. 
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Appendix A. Exact Stokes wave 
In deep water, a seventh-order Stokes expansion (for the exact equations) is 

K r] = acosO + i« 2 (l + ^a 2 + ^a 4 ) cos26> 

+ §a 3 (l + §a 2 + 3 g 2 f a 4 )cos3# + ia 4 (l + ^a 2 ) cos4# 
+ ±§§a 5 (l+^a 2 )cos50 + §a 6 cos60 + if^a 7 cos 76 + 0(a 8 ), 
g-l K %4> = a(l-ia 2 -||a 4 - 2 |§a 6 )sm6l + ia 2 (l + ^a 2 + ||a 4 )sin26l 
+ |a 3 (l + ffia 2 + f§§a 4 )sin30 + ^(l + ffa 2 ) sin40 
+ ±§§a 5 (l+ 3 ^a 2 )sin50 + fa 6 sin 60 + £a 7 sin 79 + 0(a 8 ), 

707 , 
384 ' 



l l 



g-^K^c = 1 + W + W+ 7 ||a 6 + 0(a 8 ), 



where = K{x—ct). Note that, to the leading order, the n-th Fourier coefficient is 2 1 n n n 2 a n / (n- 
1)! (this is also true for all n > 7). In the bulk of the fluid, the velocity potential is 



g-hJcp = a (l-|« 2 -n« 4 -w« 6 ) eKVsin6 + ^« 4 (l + T« 2 ) e2,ti/sin ^ 
+ ^a 5 (l+ 1 2 f a 2 )e 3f ^sin3# + ^a 6 e 4Ky sin40 

+ m a '' ' e5K% ' sin5 ° + °( a8 )' 

meaning that harmonics appear at the fourth-order only, thus justifying the ansatz (4.1) 
Note that, to the leading order, the n-th Fourier coefficient is a n+2 /n!(n — 1) for all n > 1. 



Appendix B. Cubic Zakharov's equations 

Satisfying exactly the Laplace equation and the bottom impermeability, the gravity waves 
variational formulation [36] yields the Hamiltonian [60]: 

■K = \f{gr, 2 + ^V) d 2 x, V=[cp y -Vr,.V<j ) ] y=ri . (B.l) 

Introducing a Dirichlet-Neumann operator G, such that V = G(rj)<j) (Craig & Sulem 1993, 
[12]), expanding G around r/ = and neglecting the terms beyond the quartic nonlinearities, 
the Hamiltonian (B.l) becomes 

or l f 2 



J{gr] 2 + 4> 50-O(r?50)- V- (rjVfy 

+ ±d(T] 2 V 2 4>) +d(r ] d(r ] *4>)) + \ V 2 (r/ 2 *$) ] } d 2 x, (B.2) 



with the pseudo-differential operator d = (— V 2 ) 2 tanh[(— V 2 ) 2 d]. (For one horizontal di- 
mension in infinite depth t)/ = —S)(f x ), $j the Hilbert transform.) Thus, the cubic Zakharov's 
equations (CZE) are 

Tft - Q4> = -V ■ (t/V<£) - d(r]d^) + iV 2 (r/ 2 t)0) + 8(r/5(r/D^)) + ±c)(7? 2 V 2 ^), (B.3) 
4>t + 9V = - ^(V0) 2 - (r^)V 2 - (#)^#). (B.4) 



20 



D. CLAMOND AND D. DUTYKH 



For progressive (27r/K)-periodic solutions in infinite depth, a seventh-order Stokes expansion 

is 

K r) = acos9 + ia 2 (l + §a 2 + ^a 4 ) cos26> 

+ | a 3( 1 + 41 a 2 + ^ a 4) cos3 ^ + 7_ a 4( 1+ 263 a 2-) cosAe 

+ ^a 5 (l + ^a 2 )cos5(9- ^a 6 cos 69 - ^gggg; a 7 cos 7 6 + 0(a 8 ), 

g-hjj> = a (l_l a 2_|l a 4_465 a 6) sin0 + l a 2( 1+ 3 a 2 + ^ a 4) sin20 

+ t« 3 ( 1 + !« 2 + W« 4 ) sin3 # + ^« 4 (l + ^a 2 )sin40 

+ |7a 5 (l + ^£fa 2 )sin50- J^aSin 60 - ±^a 7 sin 70 + 0(a 8 ), 

5 -5 K 3 C = 1 + i a 2 + |ia 4 + §fa 6 + 0(a 8 ), 

where the incorrect (compared to the exact expansion) coefficients and signs are displayed in 
bold face. Thus, the CZE match the exact Stokes wave up to the third-order only. Truncating 
the Hamiltonian at the order N + 1 in nonlinearities, the corresponding Stokes double series 
is correct up to the order N in the expansion parameter. None of these approximations have 
the exact asymptotic behavior (4.7) for their Fourier coefficients. 
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